Griffiths phase of the Kondo insulator fixed point 
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Heavy fermion compounds have long been identified as systems which are extremely sensitive to the 
presence of impurities and other imperfections. In recent years, both experimental and theoretical 
work has demonstrated that such disorder can lead to unusual, non-Fermi liquid behavior for most 
physical quantities. In this paper, we show that this anomalous sensitivity to disorder, as well as 
the resulting Griffiths phase behavior, directly follow from the proximity of metallic heavy fermion 
systems to the Kondo insulator fixed point. 
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Much of the recent interest in strongly correlated elec- 
tronic systems, such as heavy fermion compounds, has 
been sparked by a gamut of unusual properties found 
in experiments. Most notably, many such materials dis- 
play unusual, non-Fermi liquid (NFL) behavior which is 
often consistently observed in both the thermodynamic 
and transport measurements. While the precise origin 
of these anomalies remains a highly controversial issue, 
both the experimental and the theoretical advances have 
provided convincing evidence [jlj that disorder may be at 
the origin of such behavior, at least for certain classes 
of materials. The purpose of this paper is to clarify the 
physical content behind one of the proposed mechanisms 
for such disorder-driven NFL behavior. 

Based on an initial success in explaining the anoma- 
lous behavior of UCu5_ x Pd x by means of the so-called 
Kondo disorder model (KDM) ^|J3) , we have recently ex- 
tended the analysis to describe Anderson localization ef- 
fects in a strongly correlated environment QJ^j . It is our 
intention to show here how the emergence of NFL prop- 
erties can be described in terms of a quantum Griffiths 
phase induced by Anderson localization effects. Further- 
more, the onset of anomalous behavior occurs already 
at relatively weak disorder, an effect we ascribe to the 
proximity to the Kondo insulator fixed point. In a clean 
compound, unitary Kondo scatterers act coherently to 
create a hybridization gap and the state is the familiar 
Kondo insulator. A small deviation from unitarity leads 
to the formation of heavy fermionic quasi-particles. In 
a disordered system, however, large spatial fluctuations 
induce the appearance of random unitary scatterers or 
"Kondo insulator droplets" which are responsible for a 
strong renormalized effective disorder. These droplets, 
in turn, are regions of depleted density of states (DOS) 
which fail to quench nearby localized moments. The lat- 
ter ultimately lead to the NFL behavior. 

Heavy fermion non- Fermi liquids are characterized by 
logarithmic or weak power law divergence of the mag- 
netic susceptibility x(T) and the specific heat coeffi- 
cient j(T) — Cy(T)/T and a resistivity that behaves as 
p(T) = p + AT a , with a < 2. Deviations from normal 
Fermi-liquid behavior have also been observed in optical 



conductivity |(J, magneto- resistance |@L dynamic mag- 
netic susceptibility pL NMR and ^SR |2] measurements. 

Several mechanisms have recently been proposed to 
account for such anomalous behavior. One of them is 
provided by the proximity to an ordering transition at 
zero temperature, a quantum critical point (QCP) [||, 
where critical modes mediate unscreened long-range in- 
teractions between carriers. Indeed, magnetic ordering 
is fairly common in heavy-fermion physics and in many 
of the cleaner compounds a strong case can be made in 
favor of a QCP interpretation pf]| . However, our current 
treatment of QCP's has not produced a unified picture 
able to account for all the observed behavior of these 
systems | Q . 

Another possible route relies on the local dynamical 
frustration of exotic impurity models, where the local 
moment cannot decide with which conduction electron 
channel it will Kondo bind |l2] ]. Although a complete 
understanding of the single impurity case is available, 
there remains the question of the relevance of inter-site 
correlations jl3| . This is important because most of the 
studied systems are not in the dilute limit |l4|] . 

More recently, the idea that NFL behavior can occur 
due to the presence of disorder in a strongly correlated 
environment has been proposed. The pioneering work re- 
lied heavily on the Cu NMR line-widths of UCu5_ x Pd x 
and its temperature dependence (2). The KDM, a phe- 
nomenological model of a system with a broad distribu- 
tion of Kondo temperatures Tk, was then proposed to 
account for the data. The presence of a wide range of 
energy scales led to a picture where different spins were 
quenched at widely different temperatures. Thus, the 
singular behavior could be attributed to a few rare spins 
which remained unquenched (and therefore, highly po- 
larizable) at the lowest temperatures. These ideas were 
put on a firmer foundation by means of a dynamical mean 
field theory (DMFT) treatment of strong correlations 
and disorder, where the KDM found its natural setting 
S. The KDM had considerable success in describing 
a whole series of experiments, including thermodynamic 
responses neutron scattering jH ||, optical conduc- 

tivity and magneto-resistance |Q7 More recently, we 
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have introduced a microscopic model which incorporates 
Anderson localization effects |4||| , a feature absent from 
the KDM. We will discuss some of these recent results 
below. 

Alternatively, it has been proposed that the physics 
of these compounds is intimately tied to disorder effects 
in the proximity to magnetic ordering Jl6| . The effect of 
spatial inhomogeneities together with the tendency of the 
Kondo effect to destroy magnetic order would then lead 
to the formation of large clusters of magnetically ordered 
material close to the phase boundary but still within the 
disordered phase. Similarly, these large ordered clusters 
would be responsible for the NFL behavior. 

In all proposed disorder-based mechanisms, the low 
temperature anomalies result from a broad distribution 
of local energy scales in the system. In addition, these 
energy scales are viewed as the appropriate Kondo tem- 
peratures describing either individual spins or clusters of 
spins embedded in a metal. In all the theories, these 
Kondo temperatures assume a strong, exponential de- 
pendence on the parameters describing the fluctuator in 
question, hence their broad distribution even for moder- 
ate disorder. In the simplest KDM, the emergence of low 
Tk sites is a result of the randomness in the immediate 
environment of a given spin. In the "spin cluster" pic- 
ture, the "cluster Kondo temperature" is exponentially 
dependent on the cluster size, a quantity which can be 
expected to be large close to any magnetic quantum crit- 
ical point. Finally, in very recent work, the depression of 
Tk was proposed to be a result of Anderson localization 
effects, which lead to a reduction of the local density of 
conduction electron states required for Kondo screening. 

From a general point of view, all three proposed scenar- 
ios may be relevant and can be expected to contribute. 
A more immediate question is the relative practical sig- 
nificance of these processes and their general robustness 
with respect to different material characteristics. In this 
respect, recent experiments have provided evidence that 
the local disorder inherent to the KDM picture may not 
be sufficient to account for the observed anomalies [p"7| . 
Similarly, the emergence of large magnetically ordered 
clusters can be anticipated only in the very close vicinity 
of magnetic transitions. More importantly, it is difficult 
to imagine how the cluster picture could even come close 
to providing enough residual entropy to account for the 
observed specific heat anomalies in the physical temper- 
ature range. 

In contrast, we will show that the localization-based 
scenario provides a very robust and quantitatively rele- 
vant mechanism for the NFL behavior. Since the cor- 
responding fluctuations must be present in any moder- 
ately disordered system, this route should be of direct 
relevance to most disordered heavy fermion compounds, 
irrespective of the proximity to any magnetic ordering. 
Our key point is that the corresponding Griffiths phase 
is a direct consequence of the proximity not to any mag- 



netic, but rather to the Kondo insulator fixed point. The 
deviation from the Kondo insulator provides an energy 
scale which is universally small for any heavy fermion 
metal, since it is defined by the underlying Kondo energy. 
This observation also provides a simple explanation for 
the notorious sensitivity of HF systems to disorder, a fea- 
ture which we believe is at the origin of all the observed 
anomalies. 

We start from a disordered infinite-U Anderson lattice 
Hamiltonian 
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where, Ci a destroys a conduction electron at site i and 
spin a from a broad uncorrelated band with hopping t%j 
and fj a destroys an f-electron at site j with spin a . Since 

U — > 00, the constraint = fjafja < 1 is assumed. 

Typically, alloying introduces substitutions either in 
the f-shell sub-lattice ("Kondo holes"), e. g. Ui_ x Y x Pd3, 
or the non-f-shell subsystem ("Ligand disorder"), e. g. 
UCu5_ x Pd x . It is reasonable to assume that the former 
case should be modeled by a distribution of Efj, whereas 
the latter is expected to introduce randomness both in 
Vj and Sj. Most of our results remain unchanged irre- 
spective of the kind of disorder. 

We work within the framework of the recently intro- 
duced statistical dynamical mean field theory (SDMFT) 
It is a natural generalization of the DMFT, which 
retains the latter's treatment of local correlations while 
going beyond it by incorporating Anderson localization 
effects. It is most easily implemented on a Bethe lattice 
of coordination z (we have used z = 3 in our simula- 
tions). Each lattice site j defines an effective local action 
for the f-orbital which should be thought of as resulting 
from integrating out all the other electronic degrees of 
freedom. It is written as (U — > 00) 

5 eff = / dr fL(r) (d T + E fj ) / i)ff (r) + 

f dr f drX/l>) A ^-T')/,>(r'); (2) 
Jo Jo 
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Here, G?£t (to) is the local c-electron Green's function on 
the nearest neighbor site k with site j removed, in other 
words, a "cavity" has been created where there once was 
site j. Note how, as in the DMFT Jig ], we neglect higher 
order Green's functions in the process. G^V), on the 
other hand, is determined recursively by relating it to a 
similar quantity on the next nearest neighbor site / 
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long tails of the distribution functions. 
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Finally, the self-consistency loop is closed by requiring 
that the local self-energy Sjfe(w) be obtained from the 
solution of the effective action S^g , see Eq. (^) 0,H- 
We note that the self-consistent set of stochastic Eqs. (0- 
^) reduces to the DMFT when the limit z — > oo is taken 
(with the appropriate rescaling ~ t/y/z), in which 
case the disorder is treated on the CPA level and, there- 
fore, shows no localization effects p9| . This was a se- 
vere limitation of the previous treatment [gj which is here 
remedied. On the other hand, the non-interacting limit 
(Z7 = 0) reduces to the self-consistent theory of localiza- 
tion of Abou-Chacra, Anderson and Thouless j20|, which 
is known to exhibit an Anderson metal-insulator transi- 
tion (MIT) for z > 3. 




FIG. 1. The conduction electron wave function 9 c (r) has 
strong amplitude fluctuations due to disorder and correlate 
several local moments within a correlation length's distance. 

We would like to stress that the hybridization function 
Eq. (J3j> "seen" by each f-orbital has strong spatial fluctu- 
ations, reflecting the disorder inside a correlation volume 
enclosing several lattice sites in its neighborhood. The 
spectral information is carried by the extended conduc- 
tion electron wave function, which acts to correlate the 
different Anderson impurity problems defined by Eq. (j^) 
(see Fig. Q). This leads to a distribution of Kondo tem- 
peratures. On the other hand, the effect of this ensemble 
of single-impurity problems is to create a renormalized 
effective disorder "seen" by the conduction electrons, cf. 
Eqs. (^,||). The net effect on the relevant distribution 
functions of this highly non-local self-consistency turns 
out to be robust and universal. 

The stochastic equations (||-||) were solved by standard 
sampling techniques [po[ and provided us with the statis- 
tical distributions of the most important physical quan- 
tities. The single-impurity problem of (|[[|) was solved 
in the large-N mean-field approximation at T = 
which has the desirable feature of correctly reproducing 
the exponential form of the Kondo temperature. We have 
made sure the statistics and numerical procedures were 
accurate enough to obtain Kondo temperatures spanning 
many orders of magnitude (~ 15) in order to probe the 
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FIG. 2. Distribution of Tk showing the emergence of NFL 
behavior. Here, e,'s are distributed uniformly with width W 
and we have used z — 3, Ef = —1, V = 0.5 and fi — —0.1. 
Inset: the exponent a of Eq. (^|). The dashed line indicates 
the marginal case a = 1. 

One of our main results is the identification of a NFL 
region at relatively weak disorder Jj| . This is triggered by 
the appearance of a fraction of localized moments with 
very low Trt's, as can be seen in Fig. ||. Indeed, the Tk 
distribution exhibits power law behavior as Tk — > 



P{T K )~T\r l \ (6) 

where the exponent a varies continuously with the 
amount of disorder (see the inset of Fig. ||). This can 
be shown to lead to singular thermodynamic responses 



X (T) ~ 7 (T) 
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as has been observed in several heavy fermion alloys [E2| . 
The "marginal" case a = 1 leads to a logarithmic di- 
vergence of the same quantities. The occurrence of such 
NFL behavior in a system with a wide distribution of 
Tk's had been proposed in the context of the KDM ||. 
However, in the KDM the presence of these spins could 
only be obtained through a finely tuned choice of the bare 
disorder distribution. By contrast, in our current treat- 
ment this behavior is an unavoidable consequence of the 
spatial fluctuations of the conduction electron wave func- 
tion amplitude. Due to the extended nature of this wave 
function and the consequent correlation between Tk val- 
ues on different sites (Fig. [l]) we should expect a high 
degree of robustness and universality in these distribu- 
tions. This is indeed what we have found for different 
types of disorder distributions. We also note that the 
case depicted in Fig. |^ corresponds to conduction band 
disorder only, for which the KDM would predict no Tk 
fluctuations. 
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Diverging thermodynamic responses with disorder- 
dependent exponents due to rare regions with very low 
Tjf's are characteristic of Griffiths phases Since 
Anderson localization effects are the driving mechanism 
here, we associate this Griffiths phase to the proxim- 
ity to a disorder-driven metal- insulator transition (MIT). 
This can be checked through an examination of the typ- 
ical density of states (DOS) of the conduction electrons 
ptyp = exp{< lnpj >}; p 3 = (l/n)ImG C j(u) = 0), a 
quantity known to vanish at the MIT. We show our re- 
sults in Fig. H The Griffiths phase is observed for rel- 
atively small amounts of disorder (W/t w 1), whereas 
the MIT occurs at much higher values (W/t w 12). Sur- 
prisingly, however, the typical DOS is a non-monotonic 
function of the disorder strength (for —0.2 < p, < 0.3), 
in sharp contrast to the non- interacting case (Fig. ||). 
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FIG. 3. Localization properties of the conduction electrons 
as monitored by their typical DOS as a function of disor- 
der, for several values of the chemical potential p. The 
other parameters are the same as in Fig. ^. p = 0.1 cor- 
responds to a Kondo insulator in the clean case. While the 
non-interacting case is monotonically decreasing, the interac- 
tion-induced renormalized disorder leads to a dip for small 
values of W/t. We attribute this dip to the proximity to the 
Kondo insulator fixed point. 

In order to gain more insight into this rather non- 
intuitive behavior, we look at the effective renormalized 
disorder "seen" by the conduction electrons. It is clear 
from Eq. ([|) that, in addition to the bare disorder in the 
Efe's, scattering from the strongly correlated f-sites adds 
an effective disorder described by the quantity <5fc(u;) of 
Eq. (a). In particular, unitary scattering (S = vr/2) cor- 
responds to $(0) — > co. A diverging $(0) in the clean 
case leads to the formation of a hybridization gap and 
Kondo insulating behavior. It is precisely the appear- 
ance of the first unitary scatterers (USC's) once disorder 
is introduced which is responsible for the sharp drop in 
Ptyp seen in Fig. ^. This can be clearly seen from the 



distribution of $ fc 1 (0) as a function of disorder in Fig. |J. 
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FIG. 4. Distribution of $)7 (0) as a function of disor- 
der. The inset shows the concentration of unitary scatterers 
(<E > (7 1 (0) = 0). Same parameters as in Fig. ^, but with p — 0. 

Therefore, the following picture emerges from our re- 
sults. The clean system has a small energy scale, the 
Kondo temperature Tk, which sets the distance from 
the Kondo resonance to the Fermi energy. When disor- 
der is introduced, the positions of the Kondo resonances 
start to fluctuate and a distribution of $fc(0) ensues, giv- 
ing rise to a renormalized effective disorder. The latter 
initially broadens monotonically and USC's are quickly 
formed, having a dramatic effect on the conducting prop- 
erties. The extreme sensitivity of the system to even 
weak disorder is therefore a result of the existence of the 
small energy scale Tk, which sets the proximity to the 
Kondo insulator. Thus, the deeper in the Kondo limit, 
the more sensitive the system will be to disorder. Note 
that the value of $fc(0) is not set by Tk but is rather 
a local measure of particle-hole asymmetry (see the first 
ref. of H ) . As the disorder is further increased, the DOS 
depletion of the Kondo (pseudo-)gap is washed way and 
the concentration of USC's saturates. Finally, once the 
nearby Kondo fixed point is securely destroyed, conven- 
tional localization effects set in and the system proceeds 
to an Anderson- like MIT. That the actual MIT is not 
very much affected by the renormalized disorder is evi- 
denced by the fact that we observe it to occur at the same 
point as the non-interacting system. It is the competi- 
tion between renormalized and bare disorder which, in a 
three-stage process, leads to the non-monotonic behavior 
of Fig. | 

It is worth stressing that the emergence of the Griffiths 
phase behavior already at very weak disorder directly fol- 
lows from the described disorder renormalization as in- 
duced by the proximity of an incipient Kondo insulator. 
Indeed, this mechanism leads to a very large effective dis- 



4 



order as seen by conduction electrons, which in turn pro- 
duces the localization-induced local density of states fluc- 
tuations, and the resulting broad distributions of Kondo 
temperatures. In other correlated systems, such as the 
disordered Hubbard model |Q, the Kondo gap is not 
present, and the relevant Griffiths phase proves to be re- 
stricted to the immediate vicinity of the Mott- Anderson 
transition, in dramatic contrast to what we find here. 

In conclusion, we have elucidated the mechanism of 
the Griffiths phase observed in our studies of disordered 
Anderson lattices. Within our theory, the NFL behavior 
proves to be intimately related to the physics of disor- 
dered metals close to the Kondo insulator fixed point. 
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